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Abstract 



Measurements of open charm production cross sections in deep-inelastic ep scatter- 
ing at HERA from the HI and ZEUS Collaborations are combined. Reduced cross sec- 
tions a^ d for charm production are obtained in the kinematic range of photon virtuality 
2.5 <Q 2 < 2000 GeV 2 and Bjorken scaling variable 3 • 10~ 5 < x < 5 • 10~ 2 . The com- 
bination method accounts for the correlations of the systematic uncertainties among the 
different data sets. The combined charm data together with the combined inclusive deep- 
inelastic scattering cross sections from HERA are used as input for a detailed NLO QCD 
analysis to study the influence of different heavy flavour schemes on the parton distribution 
functions. The optimal values of the charm mass as a parameter in these different schemes 
are obtained. The implications on the NLO predictions for and Z production cross 
sections at the LHC are investigated. Using the fixed flavour number scheme, the running 
mass of the charm quark is determined. 
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1 Introduction 



Measurements of open charm production in deep-inelastic electrorfl-proton scattering (DIS) at 
HERA provide important input for stringent tests of the theory of strong interactions, quantum 
chromodynamics (QCD). Previous measurements Q3-EEH have demonstrated that charm quarks 
are predominantly produced by the boson-gluon-fusion process, — > cc, which is sensitive to 
the gluon distribution in the proton. 

The mass of the charm quark, m c , provides a sufficiently high scale necessary to apply per- 
turbative QCD (pQCD). However, additional scales are involved in charm production, e.g. the 
virtuality, Q 2 , of the exchanged photon in case of DIS and the transverse momenta, pr, of 
the outgoing quarks. The presence of several hard scales complicates the QCD calculations 
for charm production. Depending on the details of the treatment of m c , Q and p T , different 
approaches in pQCD have been formulated. In this paper, the massive fixed-flavour-number- 
scheme (FFNS) Ifl9l - |23l and different implementations of the variable-flavour-number- scheme 
(VFNS) H24H3TH are considered. 

At HERA different techniques have been used to measure open charm production cross sections 
in DIS. The full reconstruction of D or D* mesons [fTll2ll446l[T0Ul2l[T5l[T8ll. the long lifetime 
of heavy flavoured hadrons [7l-l9l [T2l[T4l or their semi-leptonic decays lfT3l are exploited. In 
general, the best signal-to-background ratio of the charm samples is observed in the analysis of 
fully reconstructed D* mesons. However, the branching ratios are small and the phase space 
of charm production accessible with D* mesons is restricted considerably because all products 
from the D* meson decay have to be measured. The usage of semi-leptonic decays of charmed 
hadrons for the analysis of charm production can profit from large branching fractions and a 
better coverage in polar angle at the cost of a worse signal-to-background ratio. Fully inclusive 
analyses using lifetime information are not hampered by specific branching ratios and are in 
addition sensitive to low transverse momenta. Among the methods used it has the largest phase 
space coverage, however it yields the worst signal-to-background ratio. 

In this paper the published data of HI HMTMHIMIl and ZEUS [|ll6lfT2lfT3|| are combined. All 
publications on data sets^l are included for which the necessary information on systematic un- 
certainties needed for the combination is available and which have not been superseded. For the 
combination, the published cross sections in the restricted phase space regions of the individual 
measurements are extrapolated to the full phase space of charm production in a coherent man- 
ner by the use of FFNS calculations in next-to-leading order (NLO). This includes the coherent 
treatment of the related systematic uncertainties. 

The combination is based on the procedure described in ll32l - [34ll . The correlated systematic 
uncertainties and the normalisation of the different measurements are accounted for such that 
one consistent data set is obtained. Since different experimental techniques of charm tagging 
have been employed using different detectors and methods of kinematic reconstruction, this 
combination leads to a significant reduction of statistical and systematic uncertainties. 

The combined charm data are used together with the combined inclusive DIS cross sections Il34fl 
to perform a detailed QCD analysis using different models of charm production in DIS. The 

'in this paper 'electron' is used to denote both electron and positron if not otherwise stated. 
2 The data taken up to the year 2000 and data taken after 2002 are referred to as HERA-I and HERA-II, 
respectively. 
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role of the value for the charm quark mass which enters as a parameter in these models is 
investigated and the optimal value of the charm quark mass parameter is determined for each 
of the QCD calculations considered. The impact of this optimisation on predictions of W ± and 
Z production cross sections at the LHC is discussed. The running mass of the charm quark is 
determined using the modified minimal subtraction scheme (MS) variant |[35l of the FFNS. 

The paper is organised as follows. In section [2] the different theoretical schemes of charm 
production are briefly reviewed. The data samples used for the combination and the details of 
the combination procedure are described in section [3j The results on the combined reduced 
cross section are presented in section @J The predictions from different QCD approaches for 
charm production in DIS are compared to the measurement in section [5J The QCD analysis is 
presented in section [6] Conclusions are given in section [7J 



2 Open charm production in DIS 

In this paper, charm production via neutral-current deep-inelastic ep scattering is considered. 
In the kinematic domain addressed, where the virtuality Q 2 of the exchanged boson is small, 
Q 2 <C M§, charm production is dominated by virtual photon exchange. The cross section may 
then be written in terms of the structure functions -F 2 cc (x, Q 2 ) and F£ c (x, Q 2 ) as 



dV 5 2vra 2 (Q 2 



-([1 + (1 - y) 2 ]Ff(x, Q 2 ) - y 2 Ff(x } Q 2 )). (1) 



dxdQ 2 xQ 4 

Here x = Q 2 /2p-q is theBjorken scaling variable and y — p-q/p-l is the inelasticity with p, q 
and / denoting the 4-momenta of the proton, photon and electron, respectively, and Q 2 = —q 2 . 
The suffix cc indicates the presence of a cc pair in the final state, including all possible QCD 
production processes. The cross section d 2 a cc /dxdQ 2 is given at the Born level without QED 
and electro-weak radiative corrections, except for the running electromagnetic coupling, a(Q 2 ) . 

In this paper, the results are presented in terms of reduced cross sections, defined as follows: 

dV 5 xQ 4 



red 



dxdQ 2 27ra 2 (Q 2 ) (1 + (1 - yf 

,2 



The contribution F£ c , originating from the exchange of longitudinally polarised photons, is 
small in the kinematic range of this analysis and reaches up to a few per cent only at high 

The above definition of F£? L Jx, Q 2 ) (also denoted as F c [|29l or F Ct si [1371 ) is suited for mea- 
surements in which charm is explicitly detected. It differs from what is sometimes used in 
theoretical calculations in which F£, L Jx, Q 2 ) ll28ll29~ll38Tl is defined as the contribution to the 
inclusive F 2 r L \(x, Q 2 ) in which the virtual photon couples directly to a c or c quark. The latter 
excludes contributions from final state gluon splitting to a cc pair in events where the photon 
couples directly to a light quark, and contributions from events in which the photon is replaced 
by a gluon from a hadron-like resolved photon. As shown in table 1 of 11291 . the gluon splitting 
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contribution is expected to be small enough to allow a reasonable comparison of the experi- 
mental results to theoretical predictions using this definition. The hadron-like resolved photon 
contribution is expected to be heavily suppressed at high Q 2 , but might not be completely neg- 
ligible in the low Q 2 region. From the point of view of pQCD it appears at 0(a 3 s ) and it is 
neglected in all theoretical calculations used in this paper. 

At photon virtualities not much larger than the charm quark mass, charm production in DIS 
is described in the framework of pQCD by flavour creation through the virtual photon-gluon- 
fusion process. Since a cc pair is being produced, there is a natural lower cutoff of 2m c for the 
mass of the hadronic final state. The non-zero mass influences the kinematics and higher order 
corrections in essentially all the HERA phase space. Therefore the correct treatment of the mass 
of charm and beauty quarks is of particular importance in the QCD analysis and determination 
of parton distribution functions (PDFs) of the proton. In the following, the different approaches 
used in the treatment of the charm quark mass in pQCD calculations are discussed. 

2.1 Zero Mass Variable Flavour Number Scheme 

In the zero-mass variable-flavour-number- scheme (ZM-VFNS) 11241 the charm quark mass is set 
to zero in the computation of the matrix elements and kinematics, and a threshold is introduced 
at Q 2 ~ m 2 , below which the charm production cross section is assumed to vanish. The 
charm quark is also excluded from the parton evolution and only three light flavours are left 
active. Above this threshold, charm is treated as a massless parton in the proton, leading to the 
introduction of the charm quark distribution function of the proton. The transition from three 
to four active flavours in the parton evolution follows the BMSN prescription [26J. The lowest 
order process for charm production in this approach is the quark-parton-model like scattering 
at order zero in a s . The running of a s is calculated using three flavours (u, d, s) below the scale 
m c , and using four or five flavours (including charm and beauty) above the respective threshold 
scales. The main advantage of this scheme is that the Q 2 evolution of the charm density provides 
a resummation of terms proportional to \og(Q 2 /m 2 ) that may be large at large Q 2 . It has been 
shown lfT5l[T8l that this approach does not describe the charm production data at HERA. 

2.2 Fixed Flavour Number Scheme 

In the fixed-flavour-number-scheme (FFNS) the charm quark is treated as massive at all scales, 
and is not considered as a parton in the proton. The number of active flavours, nj, is fixed 
to three, and charm quarks are assumed to be produced only in the hard scattering process. 
Thus the leading order (LO) process for charm production is the boson-gluon-fusion process at 
0(a s ). The next-to-leading order (NLO) coefficient functions for charm production at 0(a 2 ) 
in the FFNS were calculated in |[T9l and adopted by many global QCD analysis groups 11201 - 1231 . 
providing PDFs in the FFNS. In the data analysis presented in this paper, the prediction of open 
charm production in the FFNS at NLO is used to calculate inclusive [|T9l and exclusive Il39l 
quantities. 

In the calculations |[T9ll39l the pole mass definition 11401 is used for the charm quark mass, 
and gluon splitting contributions are included. In a recent variant of the FFNS scheme (ABM 
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FFNS) [|35l . the running mass definition in the modified minimal subtraction scheme (MS) is 
used instead. This scheme has the advantage of reducing the sensitivity of the cross sections to 
higher order corrections, and improving the theoretical precision of the mass definition. 



To 0(a s ), which is relevant for the calculation of cross sections to 0(a 2 ), the MS and pole 
masses are related by [14111 



i.e. the running mass evaluated at the scale Q = m c is smaller than the pole mass. 

2.3 General Mass Variable Flavour Number Scheme 

In the general-mass variable-flavour-number-schemes (GM-VFNS) charm production is treated 
in the FFNS in the low Q 2 region, where the mass effects are largest, and in the ZM-VFNS 
approach at high Q 2 , where the effect of resummation is most noticeable. At intermediate 
scales an interpolation is made between the FFNS and the ZM-VFNS, avoiding double count- 
ing of common terms. This scheme is expected to combine the advantages of the FFNS and 
ZM-VFNS, while introducing some level of arbitrariness in the treatment of the interpolation. 
Different implementations of the GM-VFNS are available [|25] - [3TI and are used by the global 
QCD analysis groups. 

The freedom introduced by choosing an interpolation approach prevents a clear interpretation 
of the charm mass in terms of a specific renormalisation scheme. Therefore the charm mass 
appearing in the GM-VFNS will be treated in the following sections as an effective mass pa- 
rameter, M c , of the individual interpolation models. 



3 Combination of HI and ZEUS measurements 
3.1 Data samples 

The HI ll42ll and ZEUS [43J detectors were general purpose instruments which consisted of 
tracking systems surrounded by electromagnetic and hadronic calorimeters and muon detectors, 
ensuring close to An coverage of the ep interaction point. Both detectors were equipped with 
high-resolution silicon vertex detectors: the Central Silicon Tracker B41 for HI and the Micro 
Vertex Detector [45] for ZEUS. 

The data sets included in the combination are listed in tabled] and correspond to 155 different 
cross section measurements. The combination includes measurements of charm production 
performed using different tagging techniques: the reconstruction of particular decays of D- 
mesons ll4l[6] [T0l[T2l[T5][T8ll . the inclusive analysis of tracks exploiting lifetime information [fT4l 
or the reconstruction of muons from charm semi-leptonic decays [|T3l . 

The results of the inclusive lifetime analysis Ifl4l are directly taken from the original measure- 
ment in the form of a^ d . In the case of D-meson and muon measurements, the inputs to the 



m c {Q) 




Q 2 



(3) 



m c (m c ) 2 
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combination are visible cross sections cr v i S; bin defined as the D (or /x) production cross section 
in a particular p T and rj range, reported in the corresponding publications^, in bins of Q 2 and 
y or x. Where necessary, the beauty contribution to the inclusive cross sections of D meson 
production is subtracted using the estimates of the corresponding papers. The measured cross 
sections include corrections for radiation of a real photon from the incoming or outgoing lepton 
and for virtual electroweak effects using the HERACLES program ll46ll . QED corrections to the 
incoming and outgoing quarks were neglected. All D-meson cross sections are updated using 
the most recent branching ratios ROl . 



3.2 Extraction of a^ d from visible cross sections 

In the case of D-meson and muon production, cr£f d is obtained from the visible cross sections 
c"vis,bin measured in a limited phase space using a common theory. The reduced charm cross 
section at a reference (x, Q 2 ) point is extracted according to 

cc ( n 2\ _ °"red ( X ? Q 2 ) ,a\ 

°redl x 5 l * ) ~ °vis,bin t ^ ■ V+) 

^vis,bin 

The program from Riemersma et al. [19J and the program HVQDIS [39j are used to calculate, 
in NLO FFNS, the reduced cross sections cr^ h (x, Q 2 ) and the visible cross sections cr*^ sbin , 
respectively. The following parameters are used consistently in both NLO calculations and the 
corresponding variations are used to estimate the associated uncertainties on the extraction of 

a red • 



• pole mass of the charm quark m c = 1.5 ± 0.15 GeV; 

• renormalisation and factorisation scales [if = /j r = \/Q 2 + ^m 2 , varied simultane- 
ously up or down by a factor of two; 

• strong coupling constant al f ~ ? \M z ) = 0.105 ± 0.002, corresponding to a" /_5 (M^) = 

0.116 ±0.002; 

• the proton structure is described by a series of FFNS variants of the HERAPDF1.0 
set [34] at NLO, evaluated for m c = 1.5 ± 0.15 GeV and for a™ f=3 (M z ) = 0.105 ± 
0.002. For the light flavour contribution, the renormalisation and factorisation scales are 
set to n r = fif = Q, while for the heavy quark contributions the scales of /// = n r — 

\JQ 2 + 4mg are used, with raq being the mass of the charm or beauty quark. Additional 
PDF sets are evaluated, in which the scales are varied simultaneously by a factor of two 
up or down. Only the scale variation in the heavy quark contribution has a sizeable effect 
on the PDFs. The experimental, model and parameterisation uncertainties of the PDFs 
at 68% C.L. are also included in the determination of the PDF uncertainties on a^ d . 
For estimating the uncertainties of the NLO calculations |fT9ll39ll due to the respective 
choice of the scales, a s and m c , the appropriate PDF set is used. The effects of the PDF 
uncertainties are calculated according to the HERAPDF1.0 prescription [|34l . 

3 A misprint was found in table 3 of |6|: for the rows 22 and 23 the y ranges should read 0.22 — 0.10 and 
0.10 — 0.02, respectively. Another misprint was found in table 2 of [131: the Q 2 range in the last row should be 
400 : 10000 GeV 2 . 
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The cross sections a™, bin depend, in addition to the kinematics of the charm quark produc- 
tion mechanism, also on the fragmentation of the charm quark into particular hadrons. The 
charm quark fragmentation function has been measured by HI 11471 and ZEUS [fTTI using 
the production of D* mesons, with and without associated jets, in DIS and photoproduction 
(Q 2 ~ GeV 2 ). In the calculation of a^ s bin the fragmentation is performed in the 7*-p centre- 
of-mass frame, using for the fraction of the charm quark momentum carried by the charmed 
meson a fragmentation function which is controlled by a single parameter, ll48l . The pa- 
rameter relevant for charm fragmentation into D* mesons has been determined |[TTll47ll for the 
NLO FFNS calculation for three different kinematic and jet requirements, which correspond 
approximately to three different regions of the 7*-parton centre-of-mass energy squared, s. The 
values of ax, together with the corresponding ranges in s, are listed in table [2] The fragmen- 
tation is observed to become softer with increasing s, as expected from the evolution of the 
fragmentation function. The limits on the s ranges are determined with HVQDIS by applying 
the jet requirements of the individual analysis on parton level. The olk parameters and the s 

ranges were varied according to their uncertainties to evaluate the corresponding uncertainty on 

th *~ 

vis, bin' 

Since ground-state D mesons partly originate from decays of D* and other excited mesons, 
the corresponding charm fragmentation function is softer than that measured using D* mesons. 
From kinematic considerations 11491 . supported by experimental measurements [|50l . the expec- 
tation value for the fragmentation function of charm into £)°> noD * + ? £)+ and in the mix of charm 
hadrons decaying into muons, has to be reduced by m 5% with respect to that for D* mesons. 
The values of ax for the fragmentation into ground state hadrons, used for the £)°> noD * + ; JJ+ 
and n measurements, have been re-evaluated accordingly and are reported in tabled 

Transverse fragmentation is simulated assigning to charmed hadrons a transverse momentum, 
k T , with respect to the charm quark direction, according to f(k T ) = k T exp(— 2kr/ (kr))- The 
average (kr) is set to 0.35 ± 0.15 GeV IrSTI . 

The fragmentation fractions of charm quarks into specific D mesons are listed in table [3] They 
are obtained from the average of e + e~ and ep results ll52l . The semi-leptonic branching fraction 
B(c — > /i) [|40l is also given. The decay spectrum of leptons from charm decays is taken 
from It53l 

To evaluate the extrapolation uncertainty on the extracted reduced cross section, a^ d , all the 
above parameters are varied by the quoted uncertainties and each variation is considered as a 
correlated uncertainty among the measurements to which it applies. The dominant contributions 
arise from the variation of the fragmentation function (average 3 — 5%) and from the variation 
of the renormalisation and factorisation scales (average 5 — 6%, reaching 15% at lowest Q 2 ). In 
a few cases, the symmetric variation of model parameters results in an asymmetric uncertainty 
on the cross section. In such cases, the larger difference with respect to the default cross section 
is applied symmetrically as systematic uncertainty. 

3.3 Common x — Q 2 grid 

Except for the HI lifetime analysis lfT4l . the values of a^ d for individual measurements are 
determined at the 52 (x, Q 2 ) points of a common grid. The grid points are chosen such that 
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they are close to the centre-of gravity in x and Q 2 of the corresponding a v i S) bi n bins, taking 
advantage of the fact that the binnings used by the HI and ZEUS experiments are similar. Prior 
to the combination, the HI lifetime analysis measurements are transformed, when needed, to 
the common grid (x, Q 2 ) points using the NLO FFNS calculation |[P9l . The resulting scaling 
factors are always smaller than 18% and the associated uncertainties, obtained by varying the 
charm mass, the scales and the PDFs, are negligible. For all but five grid points at least two 
measurements enter into the combination. 



3.4 Combination method 

The combination of the data sets uses the x 2 minimisation method developed for the combina- 
tion of inclusive DIS cross sections [|32l[34l . The x 2 function takes into account the correlated 
systematic uncertainties for the HI and ZEUS cross section measurements. For an individual 
data set, e, the x 2 function is defined as 

2 / , \ v ^ ('"' Y.r'f 1 "' 1 '!-' 1 '-') 2 ^ 2 

Xl p , e (m, b) = J2 J7 J — 2 + • ( 5 ) 

j \®i ,e,stat fJ> ' ) T (,"i,e,uncor ^Tl ) j 

Here //' e is the measured value of a^(xi, Q 2 ) at an (x, Q 2 ) point % and 7*' e , 5 ije>sta , t and 5i, e ,uncor 
are the relative correlated systematic, relative statistical and relative uncorrelated systematic 
uncertainties, respectively. The vector m of quantities m l expresses the values of the combined 
cross section for each point % and the vector 6 of quantities bj expresses the shifts of the cor- 
related systematic uncertainty sources, j, in units of the standard deviation. Several data sets 
providing a number of measurements are represented by a total x 2 function, which is built from 
the sum of the xLp,e functions of all data sets 

Xtot = X exp,e • (6) 

e 

The combined reduced cross sections are given by the vector m obtained by the minimisation 
of Xtot w i tn respect to m and b. With the assumption that the statistical uncertainties are 
constant and that the systematic uncertainties are proportional to m\ this minimisation provides 
an almost unbiased estimator of m. 

The double differential cross section measurements, used as input for the combination, are 
available Il54l with their statistical and systematic uncertainties. The statistical uncertainties 
correspond to <^ ie)Sta ,t in equation Q. The systematic uncertainties within each measurement 
are classified as either point-to-point correlated or point-to-point uncorrelated, corresponding 
to 7*' e and 5j e uncor , respectively. Asymmetric systematic uncertainties are symmetrised before 
performing the combination. The result is found to be insensitive to the details of the symmetri- 
sation procedure. 

In the present analysis the correlated and uncorrelated systematic uncertainties are predomi- 
nantly of multiplicative nature, i.e. they change proportionally to the central values. In equa- 
tion © the multiplicative nature of these uncertainties is taken into account by multiplying the 
relative errors 7*' e and <5i ie , U ncor by the expectation m\ 
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In charm analyses the statistical uncertainty is mainly background dominated. Therefore it is 
treated as constant independent of m l . To investigate the sensitivity of the result on the treat- 
ment of the uncorrected and, in particular, statistical uncertainty, the analysis is repeated using 
an alternative x 2 definition in which only correlated uncertainties are taken as multiplicative 
while the uncorrected uncertainties are treated as constant. In a third approach the statistical 
uncertainties are assumed to be proportional to the square root of m\ The differences between 
the results obtained from these variations and the nominal result are taken into account as an 
asymmetric procedural uncertainty and are added to the total uncertainty of the combined result 
in quadrature. 

Correlations between systematic uncertainties of different measurements are accounted for. Ex- 
perimental systematic uncertainties are treated as independent between HI and ZEUS. Extrap- 
olation uncertainties due to the variation of the charm quark mass and the renormalisation and 
factorisation scales, charm fragmentation as well as branching fractions are treated as correlated. 
All reduced cross section data from HI and ZEUS are combined in one simultaneous minimi- 
sation, through which the correlated uncertainties are reduced also at (Q 2 , x) points where only 
one measurement exists. 



4 Combined charm cross sections 

The values of the combined cross section a^ d together with uncorrected, correlated, procedural 
and total uncertainties are given in table|4] In total, 155 measurements are combined to 52 cross- 
section measurements. 

The data show good consistency, with a % 2 -value per degree of freedom, n^, of xV n dof = 
62/103, indicating that the uncertainties of the individual measurements have been estimated 
conservatively. The distributions of pulls (as defined in (341) is shown in figured! No signif- 
icant tensions are observed. For data with no correlated systematic uncertainties the pulls are 
expected to follow Gaussian distributions with zero mean and unit width. Correlated systematic 
uncertainties lead to narrowed pull distributions. 

There are in total 48 sources of correlated systematic uncertainty, including global normali- 
sations, characterising the separate data sets. The shifts and the reduction of the correlated 
uncertainties are given in tabled None of these systematic sources shifts by more than 1.2 a 
of the nominal value in the averaging procedure. The influence of several correlated systematic 
uncertainties is reduced significantly in the result. For example the uncertainties from the vertex 
analyses due to the light quark background (HI) and due to the tracking (ZEUS) are reduced 
by almost a factor of two. The reductions can be traced mainly to the different charm tagging 
methods, and to the requirement that different measurements probe the same cross section at 
each (x, Q 2 ) point. In addition, for certain kinematic regions one measurement has superior 
precision and the less precise ones follow its trend through the fit. The reduction of systematic 
uncertainties propagates to the other average points, including those which are based solely on 
the less precise measurements. 

The cross section tables of the input data sets used in the analysis (see section [3]) together with 
the full information of the correlations among these cross section measurements can be found 
elsewhere [|54l . The combined reduced cross section is presented in figure |2] as a function of 
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x, in bins of Q 2 , and compared to the input HI and ZEUS data in figure [3] The combined 
data are significantly more precise than any of the individual input data sets. This is illustrated 
in figure HI where the measurements for Q 2 = 18 GeV 2 are shown. The uncertainty of the 
combined results is 10% on average and reaches 6% in the region of small x and medium Q 2 . 
This is an improvement of about a factor of 2 with respect to each of the most precise data sets 
in the combination. 



5 Comparison to theoretical predictions 

Before proceeding to the QCD analysis including these data, it is instructive to compare them 
to various QCD predictions produced by different theory groups, for which the parameters are 
listed in table [6] This comparison tests the interplay between the gluon and/or heavy flavour 
PDFs as obtained in different schemes and the charm treatment within each scheme (section 
2), as well as the related choice of the central value for the respective charm mass parameter. 
Some of the findings in this section can be cross-related to corresponding more detailed NLO 
QCD studies in section [6[ In addition, the effect of NNLO corrections is studied here. The 
full calculation of the heavy quark coefficient functions is available at 0(a 2 ) only. The 0(0%) 
corrections listed in table[6]correspond to partial resummation corrections applied in some kine- 
matic ranges of charm production. Most predictions already contain some measured charm data 
from previous publications as input (see table[6]for details). 

In figure |5] the reduced cross section a^ d is compared with predictions of the MSTW group 
in the GM-VFNS at NLO and NNLO, using the RT standard (ZHl and the RT optimised ED 
interpolation procedure of the cross section at the charm production threshold. At NLO, the 
optimised prediction tends to describe the data better than the standard one at lower Q 2 . The 
description of the data is improved in NNLO compared to NLO. 

In figure[6]fhe data are compared to the NLO predictions based on HERAPDF1 .5 [55] extracted 
in the RT standard scheme using as inputs the published HERA-I [ 34l and the preliminary 
HERA-II combined inclusive DIS data. For the central PDF set a charm quark mass parameter 
M c = 1.4 GeV is used. The uncertainty bands of the predictions reflect the full uncertainties on 
the HERAPDF1.5 set. They are dominated by the uncertainty on M c which is varied between 
1.35 GeV and 1.65 GeV 11341 . Within these uncertainties the HERAPDF1.5 predictions describe 
the data well. The central predictions are very similar to those of the MSTW group for the same 
scheme. 

In figure [7] the data are compared to the predictions in the GM-VFNS by the NNPDF and CT 
collaborations. Both the NNPDF FONLL-A [29] and FONLL-B EB predictions describe the 
data fairly well at higher Q 2 , while they fail to describe the data at lower Q 2 . The description 
of the data at lower Q 2 is improved in the FONLL-C [30] scheme. The CT predictions 112211561 
are based on the S-ACOT-% heavy quark scheme. The NLO prediction, which is very similar to 
the FONLL-A scheme, describes the data well for Q 2 > 5 GeV 2 but fails to describe the data 
at lower Q 2 . Similar to the FONNL-C case the description of the data improves significantly at 
NNLO. 

In figure [8] the data are compared to the prediction of the ABM group in FFNS at NLO and 
NNLO, based on the running-mass scheme for both the coefficient functions and the PDFs 
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(135115711 . The uncertainties on the prediction include the uncertainties on m c , which dominate at 
small Q 2 . The predictions at NLO and NNLO are very similar and describe the data well in the 
whole kinematic range of the measurement. 

In summary, the best description of the data is achieved by predictions including partial O(a^) 
corrections (MSTW NNLO and ABM NNLO). The predictions including 0{a 2 s ) terms in all 
parts of the calculation (NNPDF FONLL C, CT NNLO and ABM NLO) as well as the MSTW 
NLO optimal scheme also agree well with the data. The largest deviations are observed for 
predictions based on 0(a s ) terms only (NNPDF FONLL A and CT NLO). As investigated in 
the next section, further differences can be partially explained by the different choices for the 
value of the respective charm quark mass parameter M c . 



6 QCD analysis 

The combined HI and ZEUS inclusive ep neutral current and charged current DIS cross sec- 
tions have been used previously to determine the HERAPDF1.0 parton density functions. In 
the current paper a combined NLO QCD analysis is performed using the reduced charm cross 
section together with the combined inclusive DIS cross sections [|34l . Since the charm contri- 
bution to the inclusive DIS cross section is sizeable and reaches up to « 30% at high Q 2 , this 
combined analysis is expected to reduce the uncertainties related to charm production inherent 
in all PDF extractions. In particular, the role of the charm quark mass m c (m c ) or the charm 
quark mass parameter M c , depending on the heavy flavour scheme, is investigated within all 
schemes discussed in section |2] 

The analysis is performed with the HERAFITTER ll58l program, which is based on the NLO 
DGLAP evolution scheme ll59l as implemented in QCDNUM ll60l . The invariant mass of the 
hadronic system is restricted to W > 15 GeV, and the Bjorken scaling variable x is limited 
by the data to x < 0.65. In this kinematic range target mass corrections and higher twist 
contributions are expected to be small. In addition, the analysis is restricted to data with Q 2 > 
Qmin = 3.5 GeV 2 to assure the applicability of pQCD. The consistency of the input data sets and 
the good control of the systematic uncertainties enable the determination of the experimental 
uncertainties on the PDFs using the x 2 tolerance of A% 2 = 1. 

The following independent PDFs are chosen in the fit procedure: xu v (x), xd v (x), xg(x) and 
xU(x), xD(x), where xU(x) = xu(x), and xD(x) = xd(x) + xs(x). Compared to the HER- 
APDF1.0 analysis, a more flexible parameterisation with 13 free parameters is used. At the 
starting scale Qo of the QCD evolution, the PDFs are parametrised as follows: 

xg(x) 

xU (x) 
xD(x) 

The normalisation parameters A g , A Uv , A dv are constrained by the sum rules. The parameter 
Bjj is set to Bjy and the constraint Ajj = Ajy(l — f s ), with f s being the strangeness fraction 



= A g x B ° ■ (I - x) c * - A' g x B v ■ (I - xf % (7) 

= A Uv x B ^ ■ (1 - xf^ ■ (1 + E Uv x 2 ), (8) 

= A dv x Bd v ■ (l-xf^, (9) 

= Ajjx 3 ^ ■ (1 - xfu, (10) 

= ApX B ° ■ (1 - x) c °. (11) 
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at the starting scale, ensures the same normalisation for the u and d densities for x — > 0. The 
strangeness fraction is set to f s = 0.31, as obtained from neutrino-induced di-muon production 
[I6TTI . To ensure a positive gluon density at large x, the parameter C' g is set to 25, in accordance 
with 11281 . 

The study involves variations of the charm mass parameter down to M c = 1.2 GeV with the 
exception of the S-ACOT-% scheme for which the M c scan starts at M c = 1.01 GeV. Since the 
starting scale Q has to be smaller than M c , the fits are performed with setting Ql = 1.4 GeV 2 
and Ql = 1.0 GeV 2 , respectively. In order to keep the variation of M c independent from a Qo 
variation, this value for Q is chosen irrespectively of the actual value of M c used during the 
variation procedure. 

The renormalisation and factorisation scales are set to Q for the VFNS and for the light quark 
contribution in the FFNS and to ^Jq 2 + 4m 2 b for the contribution of a heavy quark in the FFNS. 

For the strong coupling constant the values a s (Mz) = 0.1176 [|40l and a" /_3 (Mz) = 0.105 
with rif = 3 active flavours in the proton are used for the VFNS and for the FFNS, respectively. 
The definition of ^ and the a s order of the calculation are the same as those listed for the 
respective scheme in table[6]at NLO (ACOT-full and ZM-VFNS see S-ACOT-x). 

For each heavy flavour scheme a number of PDF fits is performed with varying M c from 1.2 
GeV to 1.8 GeV. For each fit the x 2 (M c ) value is calculated and the optimal value, M° pt , of the 
charm quark mass parameter in a given scheme is subsequently determined from a parabolic fit 
of the form 

X 2 (A«=xL»+(^^) 2 (12) 

to the x 2 (M c ) values. Here Xmm ls me X 2 value at the minimum and <r(M° pt ) is the fitted 
experimental uncertainty on M° pt . The procedure of this x 2 -scan is illustrated in figure |9] for 
the standard RT scheme when fitting only the inclusive HERA-I DIS data and when fitting 
these data together with cr^ d . The inclusive NC and CC cross sections from HERA-I alone 
only weakly constrain M c ; the value of x 2 (M c ) varies only slowly with M c . Once the charm 
data are included, a clear minimum is observed, which then determines M° pt . 

The systematic uncertainties on M° pt are calculated from the following variations of the model 
assumptions: 

• the strangeness fraction is varied in the range 0.23 < f s < 0.38. In a recent publication 
the ATLAS collaboration ll62l has observed f s = 0.5. This value of f s is also tested and 
found to have only a negligible effect on the determination of M° pt . 

• the b-quark mass is varied by ±0.25 GeV around the central value of 4.75 GeV. 

• the minimum Q 2 value for data used in the fit, Q 2 ^, is varied for the inclusive data 
from Qmin = 3-5 GeV 2 to 5.0 GeV 2 . For the charm data this variation is not applied 
because it would significantly reduce the sensitivity of the analysis on M c . However, the 
full difference on the fitted value M° pt obtained by using the cuts <5 2 lin = 3.5 GeV 2 or 
Qmin = 5 GeV 2 , is then taken as symmetric uncertainty due to the variation of Qmin- 

• the parameterisation uncertainty is estimated similarly to the HERAPDF1 .0 procedure. 
To all quark density functions an additional parameter is added one-by-one such that the 
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parameterisations are changed in equation ([8]) from A ■ x B ■ (1 — x) c ■ (1 + Ex 2 ) to 
A ■ x B ■ (1 - x) c ■ (1 + Dx + Ex 2 ) and in equations (|9l fTT]) from A ■ x B ■ (1 - a;) c to either 
A ■ x B ■ (1 — a;) c ■ (1 + Dx) or A- x B ■ (1 — x) c ■ (1 + Ex 2 ). Furthermore, the starting scale 
Q is varied to Ql = 1.9 GeV 2 . The full difference on the fitted value M° pt , obtained by 
using Qq = 1.9 GeV 2 and Q 2 , = 1-4 GeV 2 is then taken as symmetric uncertainty due to 
the variation of the starting scale Q . The total parameterisation uncertainty is obtained 
taking the largest difference in M ° pt of the above variations with respect to M opt for the 
standard parameterisation. 

• the strong coupling constant ct s (M z ) is varied by ±0.002. 

For each scheme the assumptions in the fits are varied one by one and the corresponding x 2 scan 
as a function of M c is performed. The difference between M° pt obtained for the default assump- 
tions and the result of each variation is taken as the corresponding uncertainty. The dominant 
contribution arises from the variation of Q^^, while the remaining model and parameterisation 
uncertainties are small compared to the experimental error. 

6.1 Extraction of M c opt in the VFNS 

The following implementations of the GM-VFNS are considered: ACOT full ll25l as used for 
the CTEQHQ releases of PDFs; S-ACOT-x W\ as used for the latest CTEQ releases of PDFs, 
and for the FONLL-A scheme (29J used by NNPDF; the RT standard scheme (28) as used 
for the MRST and MSTW releases of PDFs, as well as the RT optimised scheme providing 
a smoother behaviour across thresholds [31]. The ZM-VFNS as implemented by the CTEQ 
group 11251 is also used for comparison. In all schemes, the onset of the heavy quark PDFs is 
controlled by the parameter M c in addition to the kinematic constraints. 

In figure [10]the x 2 -values as a function of M c obtained from PDF fits to the inclusive HERA-I 
data and the combined charm data are shown for all schemes considered. Similar minimal x 2 - 
values are observed for the different schemes, albeit at quite different values of M° pt . In table |7] 
the resulting values of M° pt are given together with the uncertainties, the corresponding total 
X 2 and the x 2 -contribution from the reduced charm cross section measurements. The ACOT- 
full scheme yields the best global x 2 > while the best partial x 2 for the charm data is obtained 
using the RT optimised scheme. The fits in the S-ACOT-% scheme result in a very low value of 
M° pt as compared to the other schemes. 

In figure[TT]the NLO VFNS predictions for a^ d based on the PDFs evaluated using M c =M° pt of 
the corresponding scheme are compared to the data. In general the data are better described than 
when using the default values for M c and the predictions of the different schemes become very 
similar for Q 2 > 5 GeV 2 . Even the ZM-VFNS, which includes mass effects only indirectly [|25l . 
yields an equally good description of the cr^ d as the GM-VFNS, although it fails to describe 
more differential distributions of Z)* ± meson production and the lowest Q 2 bin in figure [TT] for 
which the ZM-VFNS cross section prediction is zero. 
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6.2 Impact of the charm data on PDFs 



In figure [T2] the PDFs from a 13 parameter fit using the inclusive HERA-I data only are com- 
pared with the corresponding PDFs when including the combined charm data in the fit. For both 
of these fits the RT optimised VFNS is used. The total PDF uncertainties include the parame- 
terisation and model uncertainties as described in section [6] except for the uncertainties due to 
M c , which is treated as follows: in the fit based solely on the inclusive data a central value of 
M c — 1.4 GeV is used with a variation in the range 1.35 < M c < 1.65 GeV, consistent with 
the treatment for HERAPDF1.0. For the fit including the combined charm cross sections this 
parameter is set to M° pt with the corresponding uncertainties as obtained by the charm mass 
scan for the RT optimised VFNS (table [7]). 

By comparing the PDF uncertainties obtained from the analysis of the inclusive data only and 
from the combined analysis of the inclusive and charm data, the following observations can be 
made: 

• the inclusion of a^ d in the fit does not alter the central PDFs significantly; the central 
PDFs obtained with the charm data lie well within the uncertainty bands of the PDFs 
based on the inclusive data only; 

• the uncertainties of the valence quark distribution functions are almost unaffected; 

• the uncertainty on the gluon distribution function is reduced, mostly due to a reduction in 
the parameterisation uncertainty coming from the constraints that the charm data put on 
the gluon through the jg — > cc process; 

• the uncertainty on the xc distribution function is considerably reduced due to the con- 
strained range of M c ; 

• the uncertainty on the xu distribution function is correspondingly reduced because the 
inclusive data constrains the sum xll = xu + xc; 

• the uncertainty on the xd distribution function is also reduced because it is constrained to 
be equal to xu at low x; 

• the uncertainty on the xs distribution function is not reduced because it is dominated by 
the model uncertainty on the strangeness fraction f s . 

Similar conclusions hold also for the other schemes discussed in this paper. 

6.3 Measurement of the charm quark mass 

In the VFNS discussed in the previous section the charm quark mass parameter M c does not 
correspond directly to a physical mass. This is different for the FFNS. An NLO QCD analysis is 
performed in the FFNS of the ABM group 11351 to determine the MS running charm quark mass 
m c (m c ) based on the inclusive neutral and charged current HERA-I DIS data and the charm 
cross section. For this purpose the coefficient functions as implemented in OPENQCDRAD ll20l 
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l63l are used. The strong coupling constant is evolved with setting the number of active flavours 
to n/ = 3, using a" /_3 (Mz) = 0.105. The same minimisation procedure as for the VFNS 
analysis is applied and the resulting dependence of the x 2 values from the QCD fits on the charm 
quark mass m c is shown in figure [[3] The fit of the parabolic function, defined in equation (fT2l) . 
results in a value of 

m c (m c ) = 1.26 ± 0.05 cxp ± 0.03 mod ± 0.02 param ± 0.02 Qs GeV (13) 

for the running charm mass in NLO. The errors correspond to the experimental, the model, 
parameterisation and a s dependent uncertainties. The same variations of the model and param- 
eterisation assumptions are applied as for the results presented in section [67T1 and discussed in 
section |6l The data are well described by the FFNS calculations for m c (m c ) = 1.26 GeV with 
a total x 2 = 627.7 for 626 degrees of freedom. The partial contribution from the charm data is 
X 2 = 43.5 for 47 data points. The measured value of the running charm quark mass is consistent 
with the world average of m c {m c ) = 1.275 ± 0.025 GeV [40] defined at two-loop QCD, based 
on lattice calculations and measurements of time-like processes. It also compares well to recent 
analyses J35l|64l of DIS and charm data at NLO and NNLO. 

6.4 Impact of charm data on predictions for W ± and Z production at 
the LHC 

The different series of PDFs obtained from fits to the HERA data by the M c scanning procedure 
in the different VFNSs are used to calculate cross section predictions for W ± and Z production 
at the LHC at yfs = 7 TeV These predictions are calculated for each scheme using the MCFM 
program (55) interfaced to APPLGRID (661 for 1.2 < M c < 1.8 GeV in 0.1 GeV steps, except 
for S-ACOT-x for which the range 1.1 < M c < 1.4 GeV is used. 

The predicted W ± and Z production cross sections as a function of M c for the different imple- 
mentations of the VFNS are shown in figure [14] and the values for the optimal choice M opt are 
summarised in table [8] For all implementations of VFNS a similar monotonic dependence of 
the W ± and Z production cross sections on M c is observed. This can be qualitatively under- 
stood as follows. A higher charm mass leads to stronger suppression of charm near threshold 
such that more light sea quarks are required to fit the inclusive data. More gluons are also 
needed to describe the HERA charm data. The need for more light sea quarks at the initial scale 
together with the creation of more sea quarks from gluon splitting at higher scales lead to an 
enhancement of the W ± and Z cross sections at the LHC. 

There is a significant spread of about 6% between the predictions if they are considered for a 
fixed value of M c , e.g. at M c = 1.4 GeV. Similarly, the prediction typically varies by about 
6% when raising M c from 1.2 to 1.8 GeV. However, when using the M opt for each scheme the 
spread of predictions is reduced to 1.4% for W~, 1.8% for Z and to 2% for W + production. 

This indicates that a good description of the HERA charm data correlates with a very simi- 
lar flavour composition of the quark PDFs at LHC scales, almost independent of the chosen 
scheme. The uncertainty on the W ± and Z cross section predictions due to the choice of the 
charm mass can thus be considerably reduced. However, the charm mass parameter must be ad- 
justed to a different value for each scheme, consistent with the HERA data, in order to achieve 
this result. 
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7 Conclusions 



Measurements of open charm production in deep-inelastic ep scattering by the HI and ZEUS 
experiments using different charm tagging methods are combined, accounting for the system- 
atic correlations. The measurements are extrapolated to the full phase space using an NLO 
QCD calculation to obtain the reduced charm quark-pair cross sections in the region of photon 
virtualities 2.5 < Q 2 < 2000 GeV 2 . The combined data are compared to QCD predictions 
in the fixed-flavour-number- scheme and in the general-mass variable-flavour-number-scheme. 
The best description of the data in the whole kinematic range is provided by the NNLO fixed- 
flavour-number- scheme prediction of the ABM group. Some of the NLO general-mass variable- 
flavour-number- scheme predictions significantly underestimate the charm production cross sec- 
tion at low Q 2 , which is improved at NNLO. 

Using the combined charm cross sections together with the combined HERA inclusive DIS data, 
an NLO QCD analysis is performed based on different implementations of the variable-flavour- 
number-scheme. For each scheme, an optimal value of the charm mass parameter, M° pt , is 
determined. These values show a sizeable spread. All schemes are found to describe the data 
well, as long as the charm mass parameter is taken at the corresponding optimal value. The 
use of M° pt and its uncertainties in the QCD analysis significantly reduces the parton density 
uncertainties, mainly for the sea quark contributions from charm, down and up quarks. 

The QCD analysis is also performed in the fixed-flavour-number-scheme at NLO using the 
MS running mass definition. The running charm quark mass is determined as m c (m c ) = 
1.26 ± 0.05 exp . ± 0.03 mod ± 0.02 param ± 0.02^ GeV. This value agrees well with the world 
average based on lattice calculations and on measurements of time-like processes. 

The PDFs obtained from the corresponding QCD analyses using different M c are used to pre- 
dict W ± and Z production cross-sections at the LHC. A sizeable spread in the predictions is 
observed, when the charm mass parameter M c is varied between 1.2 and 1.8 GeV, or when 
different schemes are considered at fixed value of M c . The spread is significantly reduced when 
the optimal value of M c is used for each scheme. 
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Data set 


Tagging method 


Q 2 range 


N 


£ 






[GeV 2 ] 




[pb- 1 ] 


1 HI VTX m 


Inclusive track lifetime 


5 - 2000 


29 


245 


2 HI D* HERA-I 0JJ 


D*+ 


2 - 100 


17 


47 


3 HI D* HERA-II El 


D*+ 


5 - 100 


25 


348 


4 HI D* HERA-II 031 


D* + 


100 - 1000 


6 


351 


5 ZEUS D* (96-97) @ 


D* + 


1 - 200 


21 


37 


6 ZEUS £>* (98-00) J6l 


D *+ 


1.5 - 1000 


31 


82 


7 ZEUS _D° [H2l 


jjO,noD*+ 


5 - 1000 


9 


134 


8 ZEUS d+ ma 


D+ 


5 - 1000 


9 


134 


9 ZEUS fi |H3l 




20 - 10000 


8 


126 



Table 1 : Data sets used in the combination. For each data set the charm tagging method, the 
Q 2 range, the number of cross section measurements iV and the integrated luminosity C are 
given. The data set with the £)°> noD * + tagging method is based on an analysis of D° mesons not 
originating from detectable D* + decays. Charge conjugate modes are always implied. 



s range 


a K (D*) 


a K (g.s.) 


Measurement 


s < Si 


6.1 ±0.9 


4.6 ±0.7 


[47] D*, DIS, no-jet sample 


Si < S < S2 


3.3 ±0.4 


3.0 ±0.3 


[47] D*, DIS, jet sample 


s > s 2 


2.67 ±0.31 


2.19 ±0.24 


LI 1] D* jet photoproduction 



Table 2: The a K parameters used for the longitudinal fragmentation into D* mesons and in 
ground state (g.s.) charmed hadrons. The first column shows the s range in which a particular 
value of a K is used, with si = 70 ± 40 GeV 2 and s 2 = 324 GeV 2 . The variations of a K are 
given in the second and third column. The parameter s 2 is not varied, since the corresponding 
uncertainty is already covered by the a K variations. 



/(c D*+) 


0.2287 ±0.0056 


f(c -> D+) 


0.2256 ± 0.0077 


/(c->/J°< notD * + ) 


0.409 ±0.014 


B(c -)• /i) 


0.096 ±0.004 



Table 3: Charm fragmentation fractions to charmed mesons and the charm branching fraction 
to muons. 
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g 2 [GeV 2 




X 


y 


cc 


$unc 


[% 




&cor 


{% 




fiproced % 




Stot 


{% 




2.5 


0.00003 


0.824 


0.1126 


14.0 


10.9 


0.3 


17.8 


2.5 


0.00007 


0.353 


0.1068 


9.0 


9.9 


0.2 




13.4 


2.5 


0.00013 


0.190 


0.0889 


10.0 


9.1 


2.2 




13.7 


2.5 


0.00018 


0.137 


0.0907 


9.5 


8.3 


1.4 


12.7 


2.5 


0.00035 


0.071 


0.0560 


8.7 


8.2 


0.0 


11.9 


5. 


0.00007 


0.706 


0.1466 


15.6 


10.0 


0.2 




1 


8.5 


5 


0.00018 


0.274 


0.1495 


8.4 


6i 




1.1 


lOi 




5 


0.00035 


0.141 


0.1151 


7.1 


6.7 


0.6 




9i 




5 


0.00100 


0.049 


0.0803 


9.2 




8.2 


0.6 


12.4 


7 


0.00013 


0.532 


0.2142 


8.1 


8.0 


0.2 




11.4 


7 


0.00018 


0.384 


0.1909 


10.2 


8.5 


2.1 


13.4 


7 


0.00030 


0.231 


0.1689 


4.6 


6.3 


0.4 




7i 




7 


0.00050 


0.138 


0.1553 


4.5 




5.9 


0.6 




7.5 




7 


0.00080 


0.086 


0.1156 


7.2 




6.0 


0.7 




9.4 


7 


0.00160 


0.043 


0.0925 


6.4 


7.6 


0.6 




9.9 


12 


0.00022 


0.539 


0.2983 


8.4 


7.2 


0.1 




11.1 


12 


0.00032 


0.371 


0.2852 


4.7 


6.5 


0.6 




8.1 




12 


0.00050 


0.237 


0.2342 


4.3 


5.1 


01 






6.6 


12 


0.00080 


0.148 


0.1771 


3.8 


5.7 


0.1 






6.9 


12 


0.00150 


0.079 


0.1413 


51 




6i 




0.1 






8.7 


12 


0.00300 


0.040 


0.1028 


6.1 




8.0 


0.2 




10.1 


18 


0.00035 


0.508 


0.3093 


9.2 




6.5 


1.0 


11.3 


18 


0.00050 


0.356 


0.2766 


4.7 


7.0 


01 






8.4 


18 


0.00080 


0.222 


0.2637 


3.8 


4.6 


0.6 




6.1 




18 


0.00135 


0.132 


0.2009 


3.3 


5.2 


0.0 




6.2 




18 


0.00250 


0.071 


0.1576 


31 




5.7 


0.1 






6.7 


18 


0.00450 


0.040 


0.1349 


5i 




8.0 


1.4 


10.0 



Table 4: The averaged reduced cross section of charm production, a^ d , as obtained from the 
combination of HI and ZEUS measurements. The values of the cross section are presented 
together with uncorrelated (5 unc ) correlated (<5 cor ) and procedural (5 proce d) uncertainties. The 
total uncertainty (5 to t) is obtained by adding the correlated, uncorrelated and procedural errors 
in quadrature. 
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g 2 [GeV 2 ] 


X 


y 


cc 
°rcd 


8unc[%] 


<W%] 


fiproced % 




Sm[%] 


32 


0.00060 


0.527 


0.4119 


15.1 


5.7 


0.1 


16.2 


32 


0.00080 


0.395 


0.3527 


4.3 


5.3 


0.3 


6.9 


32 


0.00140 


0.226 


0.2767 


3.9 


4.2 


0.4 


5.8 


32 


0.00240 


0.132 


0.2035 


4.8 


4.9 


0.3 


6.9 


32 


0.00320 


0.099 


0.1942 


7.1 


5.6 


0.3 


9.0 


32 


0.00550 


0.058 


0.1487 


6.9 


6.0 


0.4 


9.1 


32 


0.00800 


0.040 


0.1027 


10.7 


8.3 


0.4 


13.5 


60 


0.00140 


0.424 


0.3218 


6.1 


5.4 


1.4 


8.3 


60 


0.00200 


0.296 


0.3387 


4.3 


3.7 


0.4 


5.7 


60 


0.00320 


0.185 


0.2721 


4.7 


3.9 


0.4 


6.1 


60 


0.00500 


0.119 


0.1975 


4.7 


4.9 


0.3 


6.8 


60 


0.00800 


0.074 


0.1456 


12.0 


5.2 


0.6 


13.1 


60 


0.01500 


0.040 


0.1008 


10.6 


6.4 


0.8 


12.4 


120 


0.00200 


0.593 


0.3450 


7.1 


5.2 


0.6 


8.8 


120 


0.00320 


0.371 


0.2432 


15.9 


4.0 


2.1 


16.5 


120 


0.00550 


0.216 


0.2260 


5.2 


4.5 


0.6 


6.9 


120 


0.01000 


0.119 


0.1590 


6.6 


5.4 


0.8 


8.6 


120 


0.02500 


0.047 


0.0866 


13.7 


6.8 


1.2 


15.3 


200 


0.00500 


0.395 


0.2439 


8.1 


5.7 


0.7 


9.9 


200 


0.01300 


0.152 


0.1659 


6.7 


4.8 


0.4 


8.3 


350 


0.01000 


0.346 


0.2250 


8.8 


5.0 


4.1 


10.9 


350 


0.02500 


0.138 


0.1016 


11.2 


5.8 


5.1 


13.6 


650 


0.01300 


0.494 


0.2004 


11.1 


7.2 


1.1 


13.3 


650 


0.03200 


0.201 


0.0939 


12.4 


10.6 


0.9 


16.4 


2000 


0.05000 


0.395 


0.0622 


27.7 


14.4 


1.7 


31.2 



Table 4: continued 
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source 


data sets 


name 


snirt [tjj 


Reduction factor [%] 


x 

01 




HI vertex resolution 


r\ 1 


y4 


x 

02 


1—4 


Hi CJC etnciency 


-U.J 




c 

03 


1 


HI CST efficiency 


0.0 


98 


X 

04 


1 


B multiplicity 


-0.3 


96 


c 

05 


1-9 


c longitudinal fragmentation 


-0.9 


84 


c 

06 


1, 3, 4 


photoproduction background 


0.2 


94 


67 


1 


D + multiplicity 


0.0 


99 


X 

°8 


1 


D° multiplicity 


0.0 


AA 

99 


C 

09 


1 


D s multiplicity 


0.1 


98 


c 

010 


1 


fragmentation 


0.0 


100 


c 

Oil 


1 


HI VTX model: x-reweighting 


-0.4 


95 


c 

012 


1 


HI VTX model: pT-reweighting 


0.3 


74 


C 

013 


1 


T T 1 T 7TV J 1 / \ 1_ 4." 

HI VTX model: ?7(c)-ieweighting 


-0.3 


87 


014 




HI VTX u<is-background 


u.u 


JJ 


015 


1 


TJ 1 \ /TV A. ^ 

Hi V 1 a <p or c-quarlc 


u.Z 


OA 

yo 


C 

016 


1 


HI hadronic energy scale 


-0.1 


89 


? 

017 




HI VTX F2 normalisation 


-O.z 


y / 


c 

018 


O A 

3, 4 


HI Primary vertex fit 


0.1 


AA 

99 


r 

019 


r l A 

2-4 


HI electron energy 


A C 

0.6 


*C A 

69 


c 

O20 


2-4 


HI electron polar angle 


0.3 


11 


021 


3, 4 


TT1 | . . ' /I IPT1 A TT\ 

HI luminosity (HERA-11) 


-0.9 


OA 

80 


022 


3, 4 


TT1 i. CC /TTT~"T* A TT\ 

HI trigger efficiency (HERA-II) 


-0.3 


98 


023 


3, 4 


HI fragmentation model in MC 


-0.1 


89 


c 

024 


1—1 


JdH(JJ —> Ktttt) 


A 1 
0.1 


AO 


jf 

025 


2—0 


f{c — > JJ ) 


A 1 
0. 1 


y4. 


£ 

026 


2, j 


HI efficiency using alternative MC model 


A A 
0.4 


/3 


C 

027 


2-9 


NLO, m c 


0.5 


72 


c 

028 


A 

l-y 


NLO, scale 


1 

-1,2 


DO 


c 

029 


2-9 


c transverse fragmentation 


-U.Z 


no 

/o 


c 

030 


^ a 

2-9 


JNLU, rut 


A ^ 

0.2 


AT 

97 


r 

031 


^ a 

2-9 


NLO, a s (Mz) 


A O 

-0.2 


AC 

95 


032 


2 


tti 1 * '4. / 1 r\r\o ^ r\f\r\\ 

HI luminosity (1998-2000) 


-0.1 


97 


c 

033 


2 


TT1 *. £C /TUT* A T\ 

HI trigger efficiency (HERA-I) 


-0.2 


95 


034 


2 


HI MC alternative fragmentation 


-0.1 


70 


r 

035 


9 


ZbUS fi: d/RMUON ernciency 


A 1 

-0.1 


AO 

92 


c 

036 


9 


ZEUS /i: FMUON efficiency 


0.2 


97 


c 

037 


9 


ZEUS fi: energy scale 


0.0 


85 


£ 

038 


y 


Z/liUo /i. r T calibration 


A A 
0.0 


/2 


c 

039 


A 

9 


ZEUS /i: hadronic resolution 


A £L 

0.6 


1 1 


c 

040 


a 
9 


ZEUS /i: fr resolution 


A 1 

-0.2 


9/ 


A,, 
041 


Q 

y 


ZjiZ/Uo /i. iviv^ mouei 


n 1 
v. 1 


OD 


<>42 


9 


B(c -)• /x) 


0.1 


97 


<>43 


7, 8 


ZEUS lifetime significance 


0.5 


52 


044 


7 


f(c £)°) 


0.3 


97 


045 


8 


/(c £)+) X BR(D+ -)• Ktttt) 


-0.6 


91 


046 


7-9 


ZEUS luminosity (2005) 


-0.1 


95 


0~47 


5 


ZEUS luminosity (1996-1997) 


0.4 


96 


048 


6 


ZEUS luminosity (1998-2000) 


0.3 


90 



Table 5: Sources of bin-to-bin correlated systematic uncertainties considered in the combina- 
tion. For each source the shifts in units of standard deviations a and the reduction factor of 
the uncertainty values are given. The systematic sources corresponding to the extrapolation 
uncertainties are highlighted in bold font. The second column shows the data sets (see tabled)) 
affected by each particular source. 
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Theory 


Scheme 


Ref. 


F 2(L) 


m c 


Massive 


Massless 


a s (m z ) 


Scale 


Included 








def. 


[GeV] 


(Q 2 <m 2 ) 


(Q 2 » rnl) 


(n/ = 5) 




charm data 


MSTW08 NLO 


RT standard 


EH 


F 2(L) 


1.4 (pole) 


0{al) 


OK) 


0.12108 


Q 




MSTW08 NNLO 










approx.-O(a^) 


0(al) 


0.11707 






MSTW08 NLO (opt.) 


RT optimised 


ED 








0(a s ) 


0.12108 






MSTW08 NNLO (opt.) 










approx.-O(a^) 


0{al) 


0.11707 






HERAPDF1.5 NLO 


RT standard 


EH 


F 2(L) 


1.4 (pole) 


0(a 2 s ) 


0(a s ) 


0.1176 


Q 


HERA inclusive DIS only 


NNPDF2.1 FONLL A 


FONLL A 


E3 


n.a. 


V2 


0(a s ) 


0(a s ) 


0.119 


Q 




NNPDF2. 1 FONLL B 


FONLL B 




F 2(L) 


V2 (pole) 


0{cfi 


0(a s ) 








NNPDF2.1 FONLL C 


FONLL C 




F 2(L) 


V2 (pole) 


G(a 2 s ) 


0(a 2 s ) 








CT10NLO 


S-ACOT-x 


E2 


n.a. 


1.3 


0(a s ) 


0(a s ) 


0.118 


a/Q 2 + m 2 c 




CT10 NNLO (prel.) 




EH 


TPCC 

F 2(L) 


1.3 (pole) 


0{a 2 s ) 


0(a 2 s ) 








ABKM09 NLO 


FFNS 


E3 


TPCC 

F 2(L) 


1.18 (MS) 


0{al) 




0.1135 


^Q 2 + 4ml 


for mass optimisation only 


ABKM09 NNLO 










approx.-O(a^) 











Table 6: Calculations from different theory groups as shown in figures|5]Hl The table shows the heavy flavour scheme used and the correspond- 
ing reference, the respective F 2 (l) definition (section©, the value and type of charm mass used (equation ©), the order in as of the massive 
and massless parts of the calculation, the value of a s , the renormalisation and factorisation scale, and which HERA charm data were included 
in the corresponding PDF fit. The distinction between the two possible F 2 (l) definitions is not applicable (n.a.) for 0(a s ) calculations. 



scheme 


M opt 

[GeV] 


NC,CC cc - 
"rcd + °rcd 


xV^dp 

°rcd 


RT standard 
RT optimised 
ACOT-full 
S-ACOT-x 
ZM-VFNS 


1.50 ± 0.06 exp ± 0.06 mod ± 0.01 param ± 0.003 aa 
1.38 ± 0.05 exp ± 0.03 mod ± 0.01 param ± 0.01 as 
1.52 ± 0.05 exp ± 0.12 mod ± 0.01 param ± 0.06 as 
1.15 ± 0.04 exp ± 0.01 mod ± 0.01 param ± 0.02 Qs 
1.60 ± 0.05 exp ± 0.03 mod ± 0.05 param ± 0.01 Qs 


630.7/626 
623.8/626 
607.3/626 
613.3/626 
631.7/626 


49.0/47 
45.8/47 
53.3/47 
50.3/47 
55.3/47 



Table 7: The values of the charm mass parameter M° pt as determined from the M c scans in 
different heavy flavour schemes. The uncertainties of the minimisation procedure are denoted 
as 'exp', the model and parameterisation uncertainties are represented by 'mod' and 'param', 
respectively. Also the uncertainty due to a s variation is listed. The corresponding global and 
partial x 2 are presented per degrees of freedom n do f and per number of data points n dp , respec- 
tively. 



scheme 


a z [nb] 


a w + [nb] 


a w - [nb] 


RT standard 
RT optimised 
ACOT-full 
S-ACOT-x 
ZM-VFNS 


28.91 ±0.30 
28.85 ±0.24 
29.32 ± 0.42 
29.00 ±0.22 
28.81 ±0.24 


57.04 ±0.55 
57.03 ± 0.45 
57.84 ±0.74 
57.32 ±0.42 
56.71 ±0.40 


39.94 ±0.35 
39.93 ±0.27 
40.39 ± 0.47 
39.86 ±0.24 
39.86 ±0.25 



Table 8: NLO VFNS predictions for Z/W ± cross sections at the LHC for y/s = 7 TeV using 
the MCFM program. The calculations are based on the PDF sets extracted in the corresponding 
schemes from the HERA data using M° pt for the charm quark mass parameter. The listed cross 
section uncertainties correspond to the uncertainties on M° pt only. 



32 




Figure 1: Pull distribution for the combined data samples (shaded histogram). RMS gives the 
root mean square of the distribution. The curve shows the result of a binned log-likelihood 
Gaussian fit. 
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Figure 2: Combined reduced cross sections a^ d as a function of x for fixed values of Q 2 . The 
error bars represent the total uncertainty including uncorrelated, correlated and procedural un- 
certainties added in quadrature. 
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Figure 3: Combined reduced cross sections cr^ (filled circles) as a function of x for fixed 
values of Q 2 . The error bars represent the total uncertainty including uncorrected, correlated 
and procedural uncertainties added in quadrature. For comparison, the input data are shown: 
the HI measurement based on lifetime information of inclusive track production is represented 
by closed squares; the HI measurements based on reconstruction of D* mesons in HERA-I / 
HERA-II running periods are denoted by filled up (down) triangles; the ZEUS measurement 
using semileptonic decays into muons is represented by open circles; the ZEUS measurements 
based on reconstruction of D* mesons are depicted by open squares (open triangles) for data 
collected in 1998-2000 (1996-1997) years; the ZEUS measurements based on reconstruction 
of D° (D + ) mesons are shown by open diamonds (crosses). For presentation purpose each 
individual measurement is shifted in x. 
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H1 and ZEUS 
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Figure 4: Combined reduced cross sections cr£f d (filled circles) as a function of x for Q 2 = 18 
GeV 2 . The error bars represent the total uncertainty including uncorrelated, correlated and 
procedural uncertainties added in quadrature. For comparison, the input data are shown. For 
further details see figure [3j 
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Figure 5: Combined reduced cross sections cr^ed (filled circles) as a function of x for fixed values of Q 2 . The error bars represent the total 
uncertainty including uncorrected, correlated and procedural uncertainties added in quadrature. The data are compared to MSTW predictions 
at NLO (left panel) and NNLO (right panel). The predictions obtained using the standard (optimised) parametrisation are represented by the 
shaded bands (solid lines). The uncertainties for the optimised parametrisation are not yet available. 
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Figure 6: Combined reduced cross sections a^ d (filled circles) as a function of x for fixed values 
of Q 2 . The error bars represent the total uncertainty including uncorrelated, correlated and 
procedural uncertainties added in quadrature. The data are compared to the NLO predictions 
based on HERAPDF1 .5 extracted in the RT standard scheme. The line represents the prediction 
using M c = 1.4 GeV. The uncertainty band shows the full PDF uncertainty which is dominated 
by the variation of M c . 
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H1 and ZEUS H1 and ZEUS 




Figure 7: Combined reduced cross sections er^ed (filled circles) as a function of x for fixed values of Q 2 . The error bars represent the total 
uncertainty including uncorrelated, correlated and procedural uncertainties added in quadrature. The data are compared to predictions by the 
NNPDF group (left panel) and the CTEQ group (right panel). The predictions from NNPDF2.1 in FONNL-A, -B and -C schemes are shown 
with their uncertainties (bands with different hatch styles). The CT10 NLO prediction with its uncertainties is shown by the shaded bands. The 
uncertainties on the CT10 NNLO (prel.) predictions are not yet available. 
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Figure 8: Combined reduced cross sections a^ d (filled circles) as a function of x for fixed 
values of Q 2 . The error bars represent the total uncertainty including uncorrected, correlated 
and procedural uncertainties added in quadrature. The data are compared to predictions of the 
ABM group at NLO (hashed band) and NNLO (shaded band) in FFNS using the MS definition 
for the charm quark mass. 
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Figure 9: The values of % 2 (M C ) for the PDF fit to the combined HERA DIS data in the RT 
standard scheme. The open symbols indicate the results of the fit to inclusive DIS data only. 
The results of the fit including the combined charm data are shown by filled symbols. 
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Figure 10: The values of x 2 (M c ) for the PDF fit to the combined HERA inclusive DIS and 
charm measurements. Different heavy flavour schemes are used in the fit and presented by lines 
with different styles. The values of M° pt for each scheme are indicated by the stars. 
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Figure 1 1 : Combined measurements of a^ d as a function of x for given values of Q 2 is shown by 
filled symbols. The error bars represent the total uncertainty including uncorrected, correlated 
and procedural uncertainties added in quadrature. The data are compared to the results of the fit 
using different variants of the VFNS (represented by lines of different styles) choosing Mc = 
M° pt . The cross section prediction of the ZM-VFNS vanishes for Q 2 = 2.5 GeV 2 . 
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Figure 12: Parton density functions x ■ f(x,Q 2 ) with / = g,u v ,d v ,u,d,s,c for (a) valence quarks and gluon and for (b) sea anti-quarks 
obtained from the combined QCD analysis of the inclusive DIS data and a^ d (dark shaded bands) in the RT optimised scheme as a function of 
x at Q 2 = 10 GeV 2 . For comparison the results of the QCD analysis of the inclusive DIS data only are also shown (light shaded bands). The 
gluon distribution function is scaled by a factor 0.05 and the xd distribution function is scaled by a factor 1.1 for better visibility. 
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Figure 13: The values of x 2 f° r the PDF fit to the combined HERA DIS data including charm 
measurements as a function of the running charm quark mass m c (m c ). The FFNS ABM scheme 
is used, where the charm quark mass is defined in the MS scheme. 
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Figure 14: NLO predictions for (a) W + , (b) W~ and (c) Z production cross sections at the 
LHC for y/s = 7 TeV as a function of M c used in the corresponding PDF fit. The different lines 
represent predictions for different implementations of the VFNS. The predictions obtained with 
PDFs evaluated with the M° pt values for each scheme are indicated by the stars. The horizontal 
dashed lines show the resulting spread of the predictions when choosing M c = M° pt . 
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